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We establish bounds on the decay of time-dependent multipoint correlation func¬ 
tionals of one-dimensional quasi-free fermions in terms of the decay properties of 
their two-point function. At a technical level, this is done with the help of bounds on 
certain bordered determinants and pfaffians. These bounds, which we prove, go be¬ 
yond the well-known Hadamard estimates. Our main application of these results is 
a proof of strong (exponential) dynamical localization of spin-correlation functions 
in disordered XY -spin chains. 
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1 Introduction and main results 

1.1 Quasi-free Fermions on the lattice 

Systems of (spinless) fermions are described on the fermionic Fock space over a separable 
single-particle Hilbert space T-L. Fermionic annihilation and creation operators c(/) and c*(/) 
associated with / G act on and satisfy the canonical anticommutation relations (CAR), 

i.e., for any f,geH 

{cif),c*{g)} = c{f)c* {g) + c* {g)c{f) = {f,g)l (1.1) 

and all other anticommutators vanish. By A{'H), we denote the C'*-algebra generated by the 
identity, 1 , and the operators c(/) and c* {g) for all /, 5 G this is the CAR algebra associated 
to Ti. A state uj on A{'H) is any normalized, non-negative linear functional. The state u, or 
more generally any linear functional, is called quasi-free (or: quasi-gaussian, cf. E 0) if all 
correlation functions are computed from Wick’s theorem. More precisely, ui is quasi-free if for 
any n G N and any collection fj ,gj G "H giving rise to operators 

Cj ■■= c{fj) + c*{gj) , j = 1..., 2n , (1.2) 

we have u{Ci ■ ■ ■ C2n-i) = 0 and 

u{Ci ■ ■ ■ C 2 n) = ^'sgnTT W (C'^(l)C'^(2)) ' ' ' W (C'^( 2 n-l)C' 7 r( 2 n)) [ =: pfC] • ( 1 - 3 ) 

TT 

The (primed) sum is over all permutations vr which satisfy 7r(l) < 7r(3) < ••• < 7r(2n— 1 ) and 
7r(2j — 1 ) < 7r(2j) for all 1 < j < n. The right side is known as the pfaffian of the triangular 
array 

C := {^{CjCk))i^j^j^^ 2 n 

which equivalently may be identified with a skew-symmetric matrix, = —C. 

In general, fermionic quasi-free states are characterized by uj{CjCk) = {^ ^ 

terms of a one-particle density matrix F on 0 "H. The latter satisfies 0 < F < 1 and is of 

the from F = 1 - p) ~ ~ ^ simple 

ubiquitous subclass of quasi-free states a = 0 and hence they are uniquely characterized by a 
one-particle density operator 0 < p < 1 on , 

(CjCk) = {fk, Q Qj) + {fj, (1 - q) Qk) ■ (F 5 ) 

Examples of such one-particle density operators which will play some role in this paper include: 

1 . the thermal equilibrium state corresponding to a single-particle Hamiltonian H on T-L with 
inverse temperature /3 > 0 and chemical potential /r G M: 

( 1 . 6 ) 

In the limit /3 —)• oo, this turns into the fermionic ground state, g = i.e., the 

spectral projection of H corresponding to energies strictly below g. 
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2 . those of the form 


(1.7) 


Qa — \4>j){4’j 

with a € { 0 , 1 }^ and an orthonormal basis of eigenvectors of a single-particle 

Hamiltonian H on T-L. For a with the property that aj = 1 whenever Aj < and Oj = 0 
otherwise, this state coincides with the ground state of the non-interacting Fermi system 
with Fermi energy /r. Other choices of a correspond to excited states. 

We will be interested in lattice fermions for which Tf = - and primarily with d = 

1. In the lattice case, for any state oJq of the form (11.51) . the canonical orthonormal basis of 
vectors <5^ £ {lA) localized at the lattice sites ^ £ gives rise to time-dependent multipoint 
correlation functions with determinantal structure, i.e., for xi,..., yi, ■ ■ ■ , j/n G and 

•Sl) ■ ■ ■ ) Sn ) f 1) ■ ■ ■ ) fn £ K- 

= det . (1.8) 

Here for a given single-particle Hamiltonian H on ^^(Z'^), a simple (free) dynamics is imple¬ 
mented on the corresponding CAR algebra Ad ■= Al(£^(Z'^)) through 

Tt(c(/)) := c(e-'*^/), f£M; (1.9) 

more general dynamics will be considered later. (The derivation of (11.81) from (11.51) proceeds by 
setting /i = ••• = /„ = = ... = = 0 and gj = fn+j = 

for j = 1,... n and using the fact that the arising triangular array (11.41) has a block structure in 
which case the pfaffian reduces to a determinant.) 


1.1.1 Determinant bound 

In its simplest form, the basic question of this paper can already be formulated in this setting. 
Suppose the two-point function decays, i.e., for some C, g £ (0, oo) and all Xj,yk £ Z'^ 

sup < (Z, ( 1 . 10 ) 

with some monotone increasing K : [0,oo) —>• [0,oo), what can one say about the decay of 
multipoint correlation functions such as (11.81) ? Before describing our answer to this question, let 
us first clarify some points. 

Exponential decay of the two-point function would correspond to the choice K{t) = r. 
Slower decay rates such as algebraic decay may be accommodated by another choice such as 
K{t) = ln(l -I- r). Our basic assumption throughout this paper is: 

Assumption: K : [0,oo) —)■ [0,oo) is monotone increasing and there exists some /tq € 

( 0 , oo) such that 

OD 

/(/io) := J](l + ^)e-^°^W <oo. ( 1 . 11 ) 

1=0 
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Since the multipoint correlation functions depend on configurations x = (xi ,... ,Xn) G 
and y = {yi,... ,yn) G we need to specify our notion of distance for configurations. Since 
we focus here on one-dimensional systems, d = 1, we may restrict without loss of generality to 
fermionic configurations which are naturally ordered: 


Xi < X2 < ■ ■ ■ < Xn, yi < y2 < ■ ■ ■ < Vn- ( 1 - 12 ) 

Note that there is only one particle per site, since we are dealing with (spinless) fermions. In this 
situation, the notion of distance we will adopt is: 

D{x,y) := max \xj — yj\ . (1-13) 


Some remarks are in order: 

1. Since we deal with configurations of indistinguishable particles, one might wonder whether 
this distance is invariant under relabeling of particles. This immediately follows from the 
fact that D{x, y) = min.^- maxjg{i „} \xj — yjrfj)! where the minimum is over all per¬ 
mutations TT of n elements. 

2. There are other notions of distance then (11.131) which are natural for indistinguishable 

particles. For example, Di{x,y) := minTr \^j ~ d7r(j)l option which in our 

one-dimensional situation turns out to be Di{x,y) = Ylj \^j ~ djl ordered 

fermionic configurations (11.121) . 

One result of this paper is the following: 

Theorem 1.1 (Determinant bound). Let I C M and p{s,t) be a family of uniformly bounded 
operators on £^(Z) with ||/9(s, t)\\ < 1 for all s,t £ 1. If there is some C G (0, oo), p, G {po, oo) 
such that for all x,y £ Z.' 


sup \{6x,p{s,t)dy)\ < Ce ^ 1 ), (1.14) 

s,tGl 


then for any n G N and any pair of fermionic configurations x = (xi,..., Xn),y = {yi, • • •, Pn) G 

Z".- 


sup 


det [{dxj , p{sj,tk)dyf)) i<j k<n 


< 8maK{CI{po), s/CI{po)} exp 


. (1.15) 


Some remarks apply: 

1. The theorem does not require p{s,t) to be a reduced density operator - only the norm 
bound ||/9(s,f)|| < 1 is essential. In particular, the choice p{s,t) = with 

some reduced density operator 0 < p < 1 is admissible, cf. (11.81) . 


4 











2. Exponential decay of correlations, i.e., (11.141) with K{t) = r, is known to occur for ther¬ 
mal states of single-paiticle systems. More precisely, in case p = p{H) with a self-adjoint 
operator H on and any function p : M —)• C, which has an analytic extension to the 

strip I lmz\ < rj on which it is bounded by ||p||oo> obeys for all Xj,yk G 

\{6^.,p{H)6y.)\ < 18^2 IIpIIoo (1.16) 

for any p > 0 such that b{p) := sup^/g^d |(%) “ 1) < ^/2, cf. ||2l 

Theorem 3]. 

In particular, this applies with p < 7r//3 to the Fermi distribution function, p{H) = (1 -|- 

1.1.2 Disordered case 

Part of the motivation for Theorem II.II stems from the analysis of disordered systems. Here 
the single-particle Hamiltonian H (as well as the one-particle density matrix g) is a weakly 
measurable map from some probability space (fl, S,P) into the space of self-adjoint operators 
on The most prominent examples of such H are Anderson-type opeartors, i.e., discrete 

Schrddinger operators of the form H = —A + V where the multiplication operator V is given 
by independent and identically distributed (iid) random variables (t'g) associated to ^ G The 
spectral theory of such operators is quite well studied - the main feature being the existence of 
a localized phase ifTTlIISl I^. For convenience of the reader, let us summarize some facts which 
are important in the following: 

1. A dynamical characterization of localization involves the eigenfunction correlator which 
is defined as the total variation measure associated with x, y G Z'^: 

Q{x,y;I):= sup 5y)\ (1.17) 

/eL“(M) 

ll/l|oc<l 

where Pi{H) denotes the spectral projection on / C M. Strong exponential dynamical 
localization in I then refers to the bound: 

E[Q(x,y;/)] < 67(1.18) 

for some C, p G (0, oo) and all x, y G Z*^. It implies that the spectrum of H is almost- 
surely pure point. In case the latter is simple given by Ai < A2 < • • •, the eigenfunction 
correlator is then Q{x, y; I) = Z]Ajeo-(H)n/ \4>j{x)\\4>j{y)\< given in terms of the normal¬ 
ized eigenbasis {(pj} of H. 

2. Under some reasonable assumptions on the random operator H, strong exponential dy¬ 
namical localization (11.181) is known to occur in case d = 1 throughout the spectrum, i.e. 
with / = M; for details, see ifTTl ISll and references therein. 

The assumption of the following theorem hence applies to such random operators H on (Z) 
and all (time-evolved) one-particle operators p{s,t) = p{H) which are bounded 

functions of the Hamiltonian such as, for example, thermal-states (11.61) up to /3 = 00 or eigen¬ 
states (11.71) related to H. 
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Theorem 1.2 (Strong dynamical localization I). Let /CM and consider a family of random 
operators p{s,t) on /^(Z) with ||/9(s,t)|| < 1/or all s,t £ I which exhibit localization in the 
sense that for some C, p £ (0, oo) and for all x,y £ Z; 


E 


sup I(4,p(s,t) 6y)\ 
S,t£l 




(1.19) 


Then for any n £N and any pair of configurations x = (xi,..., Xn),y = (yi,... , y-n) S 


E 


sup 


det {{5o,^,p{sj,tk)6y^)) 


l</,fc<7T. 


^ SmaxIC, \/C| / 

- (l-e-M)2 exp(-f/l(x,y)) . (1.20) 


Several remarks apply: 

1. The case p{s,t) = in Theorem 11.21 includes a statement on the determinant of 

the time evolution operator projected to a pair of configurations x,y £ Z"'. In 

10 this quantity arises in the analysis of an error-correcting code for a one-dimensional 
chain of Majorana fermions. More precisely, the random one-particle Hamiltonian H 
on ..., //}) is dubbed in 0 multipoint dynamical localized if there are constants 

C, p £ (0, oo) such that for all n < sufficiently large: 


sup E 
teK 


det {{dxj 







( 1 . 21 ) 


for all configurations x, y G Z” with i9i(x, y) > N/8. 

Eq, (11.201) is weaker in case 1 <C n <C A^, since Di{x,y) < nD{x,y). It is an inter¬ 
esting open question whether (11.211) holds in the regime of strong-dynamical one-paiticle 
localization (11.181) . 

2. The exponential decay in (11.191) . and then subsequently in (11.201) . can also be replaced 
by a slower or faster decay (captured by K) as in Theorem 11.11 It is an interesting open 
question whether the above Theorems Il.ll and fOl can be generalized to higher dimensions. 


1.2 Majorana Fermions 

For lattice fermions one may associate to each site x £ if a pair of Majorana fermions: 

a+ := c*(4) -h c(4), af :=i (c*(4) - c(4)) , (1.22) 

These operators are self-adjoint, {at)* = at, and satisfy {at)"^ = 1 for both They 

obey the anticommutation relations 

{af > 4} = 24,3/<)#,b 1- (1-23) 

The free dynamics (11.91 ) generated by a single-particle Hamiltonian H carries over to the Majo¬ 
rana fermions 

W :='Tt (af) • (1-24) 
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Given such a dynamics and a quasi-free state oo the dynamical multipoint correlation functions 
are of the form 

^ [afl (G) • • • (f2n)) = pf (w (at^ (tfc)) ) 

In the subsequent theorem, we envoke the following definition. 

Definition. A pair (cj, r) of a functional io and automorphisms r = {rtltgK on the CAR algebra 
Ad ■= is called quasi-free, if the Wick relation (11.251) holds for all n € N at all 

(xi,#i),..., (X2n,#2n) S X {±}, and all times (G, .. . ,t2n) G 

As before, our main concern will be the decay rate of such multipoint correlation functions, 
given information about the decay of the two-point function. Since the multipoint correlation 
function involves a collection x := (xi,... ,X 2 n) G of points, we again first need to 

quantify the relevant notion of distance for this collection. In the one-dimensional situation, 
d = I, the points may be ordered without loss of generality 


Xl < X2 < ■ ■ ■ < X2n-l < X2n ■ (1-26) 

Note that in contrast to (11.121) these points are not necessarily distinct since Xj € Z may carry 
two Majoranas: one with # = -I- and one with ^ We will call (xi, #i), ..., (x2ru #2n) G 
Z X {zb} a Majorana configuration if these tupels are distinct for all j G {1,..., 2n}. A natural 
notion of distance for such an ordered Majorana configuration is 

r(x) := max |x2j — X2j-i| • (1-27) 

In this context our first main result then reads as follows: 


Theorem 1.3 (Pfaffian bound). Let be a quasi-free pair on the CAR algebra Ai and 

assume that ui is a bounded functional, i.e., there is some Mq G (0, oo) such that |a;(A)| < 
Afoll A|| for all A G Ai. Let / C M and suppose there is some C G (0, oo), p G {po, oo) such 
that for all x,y G Z.' 


max sup 
#,be{±} s,te/ 


ui 





(1.28) 


Then there is some C = Cfpfj such that for any n G N, and any Majorana configuration 
(3:1, ■ ■ ■ 1 {X2ni i^2n) G Z X {zb}.' 


sup 


pf (^cj (^a*J 


l<j<k<2n 


<Mo(Z'(/io)e-(«)^(’'(^)). 


(1.29) 


Let us stress that we do not assume here that w is a state: it neither needs to be non-negative 
nor normalized. Only its Gaussian nature and boundedness are essential. 

Similar to before, this bounds carries over to the random case. 
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Theorem 1.4 (Strong dynamical localization II). Let (a;, r) be a random quasi-free pair on 
the CAR algebra Ai and assume that oj is a bounded functional, i.e., there is some random 
Mq E (0, oo) such that |cc(74)| < Mo||A||/or all A E .Ai- Let / C M and suppose there is some 
(non-random) C, p, £ (0, oo) such that for all x, y E Z; 


max E 
#,be{±} 


sup 

s,tGl 


U) (of (t)a^(s)) 




(1.30) 


Then there is some C' = C'{p) E (0, oo) such that for any n E N, and any Majorana configu¬ 
ration (xi, #i), . . . , {X2n, #2n) E Z X {zb}.' 


E 


TT 

i6/2n Mq 


pf (cj [at-{tj)a 


ilk 

■j 


l<j<k<2n 


< C'{p) 


(1.31) 


Let us conclude with two remarks: 


1. It is straightforward to see from the subsequent proof that we may extend Theorem 11.41 
to the case in which uj depends on an additional parameter a. If one assumes uniform 
exponential decay in the sense that (11.301) holds with an additional supremum over a inside 
the expectation, then (11.311) holds with an additional supremum inside the expectation. 


2. Again, the above theorem has a straightforward generalization to the case that the two- 
point function decays at a rate given by K. 


2 Time-dependent correlations in random XY spin chains 


Our main application of Theorem 11.41 concerns the correlation functions of a (random) spin- 
i chain. More precisely, we consider an anisotropic spin chain of length E N with the 
Hamiltonian 

N-l N 

Sn = -Y^ /r^[(l + + (1 - ~ (2.1) 

S=i ?=i 

which acts on the Hilbert space TLn = The real-valued sequences{/r^}, { 7 ^}, and 

{i/g} are the parameters of the model which can be physically interpreted as an interaction 
strength, the anisotropy, and an external magnetic field in 3-direction, respectively, and 


a 


1 



0 —z \ . 3 

. ^ j, and .T 



( 2 . 2 ) 


denote the Pauli matrices. By the subscripts ^ E {!,... , A^}, we embed these matrices into 
B{TLn), i-e-, o-f = 1 ( 8 ) • • • ( 8 ) 11 < 8 ) u"' < 8 ) 11 < 8 ) • • • ( 8 ) 11 for any w E {1, 2 ,3} with u"' appearing in 
the ^th factor. 

The dynamics generated by the Hamiltonian Stv is the one-parameter group of automorphisms 
on given by 

rf (A) = Ae"***^^ for all A E tv) and f E M . 


(2.3) 














We are interested in dynamic correlations between general single-site observables. More con¬ 
cretely, for any 1 < ^ < denote by the set of observables with support {^}. With 
1 < ^ < ry < fixed, let A G and B G ■ We consider 

{Tt{A)B) - {Tt{A)){B ), with (•) := tr {p{Sn) (•)) • (2-4) 


The states (•) are described in terms of their density matrices p{Sn) > 0. We will mainly 
consider either eigenstates or thermal states associated to Sn, i-S-, 


p{Sn) 


|4'a)('I'al ) eigenstate of Sn with label a, 

^-PSn j q-PSn ^ thermal state with inverse temperature (3. 


(2.5) 


In order to distinguish the two cases, we will sometimes include a subscript a (in case of an 
eigenstate to be described below) or /3 (in case of a thermal state). Note that since these are 
expectations in a state whose density matrix commutes with Sat, it is clear that they are time 
invariant, i.e. {Tt{A)) = (^). To calculate the correlations (12.41) . we first expand the single-site 
observables in terms of a basis. Any A G -4{^} can be written as: A = aol + aia^ + a 2 (T| -|- 
= Ylw=o and we have set cr® = 1 for convenience. As a result, 

{r,{A)B} - {A)(B) = y; o„6„, - (»;”)«')) . (2.6) 

w,w'=l 


In order to estimate these correlation functions, we relate them to correlations of free Majorana 
fermions using the well known Jordan-Wigner transformation ifTSlfTTl . 


2.1 Jordan-Wigner transformation in terms of Majorana Fermions 


The operators 


= (t| and for all 2 < .^ < A^ , 

a~[ = —af and = —af ■ ■ ■ for all 2 < ^ < A^ , (2.7) 

are self-adjoint, (a|^)* = a^, and satisfy (a|^)^ = 1 as well as the anti-commutation rules (11.231) 
for Majorana fermions. A short calculation also shows that ia '^= cr| and that the Hamilto¬ 
nian coincides with the following quadratic form 

Sn = Ia^HnA ( 2 . 8 ) 


in terms of the vector A = {af, ,..., o^)^. The 2N x 2N coefficient matrix Hn is 

self-adjoint and of Jacobi block-form (with blocks composed of Pauli matrices): 


( -/ri5(7i) ^ 

-pN-\S{^N-\) 

\ -pN-lS{-iN-\)* VNa'^ / 


(2.9) 
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where <S'( 7 ) := The operator acting on ^^({1,... , A^}; C^) will he referred to 

as the single-particle Hamiltonian. Let us briefly comments on some of its properties: 

1. Since Hqq =: iKqq with Kjq real and skew symmetric, the spectrum of Hjq is symmetric 
about the origin, i.e., = {±Ai, ±A 2 ,..., ±AAr} with 0 <Ai<A 2 <---<AAr 

denoting its non-negative eigenvalues. 


2. The unitary transformation u ^ rotates the spin matrices, u*a^u = and 

u*a‘^u = a^. Through U := 0^i u one may lift this rotation to a local transformation 

on ^^({1,... , A^}; C^) ~ 0jLi Under this transformation, the Hamiltonian (12.91) 
turns into a Jacobi matrix block matrix U*Hj\fU, in which the variables are on the 
diagonal. 

Performing another change of variables, in which we permute the indices in U*HnU, the 
Hamiltonian is seen to be unitarily equivalent (denoted here by ~) to the block matrix 


Hn^ 




with A 


m ■■■ 

HN-l 

\ fJ-N-1 —I'N / 


( 2 . 10 ) 


and 


/ 0 71^1 


\ 


B = 




( 2 . 11 ) 


\ —'iN-lfJ-N-l 0 J 

In the isotropic case, i.e. 7 ^ = 0, and if the spin coupling is homogeneous, i.e. /rg = /r for 
all the Hamiltonian thus reduces to (two copies of) a discrete Schrddinger operator on 
^^({ 1 ,..., A^}) with hopping /r and potential given by 


To diagonalize Sat, we make a Bogoliubov transformation. More precisely, let O be the real 
orthogonal 2N x 2N matrix which brings the skew-symmetric matrix Kqq into its canonical 
block form, 

(9iTArO^ = A = Aj , where Aj := f ^o)' (2.12) 

j=i J / 

Regarding this as a change of variables and recalling the quadratic form relation (12.81) . it is 
natural to define 


B := OA and then label H = ( 6 )^, 6 ^ , • • • , 6 )^, b^y)^ , (2.13) 

in analogy to A. The Hamiltonian is then in its canonical form in terms of these 6 -operators: 

1 ^ 

57V = - A'^HnA = Xj ib+bj . (2.14) 

i=i 

Let us summarize some basic properties of these operators: 
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1. Since O is a real orthogonal matrix, the algebra of Majorana fermions is preserved under 

this transformation, i.e., {b^)* = bf, = 1 and = 26j^kb#,\, for both 

#, b e {±} and all N}. 

2. For j G ,A^}, the operators ib^bj are self-adjoint and pairwise commute. Since 

{'ib~jb~)‘^ = 1 their eigenvalues are ±1. They measure the individual fermion parity. 
More precisely, the Fermi creation and annihilation operators 

\ ^ 

corresponding to these Majorana modes, satisfy 

-t = ib+bj . (2.16) 

The Hamitonian Sn commutes, [Pn, S'at] = 0, with the total fermion parity 

Pn ■■= ib^b];^-■-ibfb^ . (2.17) 

Since the latter is self-adjoint P^ = Pn and satisfies = 1, it also has eigenvalues ±1. 
The orthogonal transformation (12.131) preserves the fermion parity operator, i.e. 

Pn = det O ■ ■ ■ ■ iafa^ = det O ■ a% .. .af . (2.18) 

(This follows most easily by restricting wlog to the case det O = 1, for which the or¬ 
thogonal transformation can be implemented on the Hilbert space by a unitary dynamics 
generated by a quadratic Hamiltonian which commutes with Pn, cf. |[8l.) 

Since the spin Hamiltonian is quadratic and diagonal in ib'^bj, a number of important con¬ 
sequences follow: 

1. The spectrum of Sn can be completely described in terms of the joint eigenstates of the 
collection of the operators ib^bj. To do so, we start from the unique normalized vector 
bl ^ T-Ln defined by 'ipjbl = 0 for all 1 < j < Next, for a = (ai, 02 , • • • , axr) G 
{0,1}^, the vectors 

= (2.19) 

form an orthonormal basis of "Htv- In fact, they are also eigenvectors of Sat: 

Sn^o. = Ea^a with Ea = 2 ^ Aj - P (2.20) 

j-.aj=l 

where E = Aj stands for the negative ground-state energy. The fermion parity of 

these eigenstates is 

PN'b^ = (-l)^.-i (2.21) 

for all a G {0,1}^. The ground-state is unique and given by 'I'(o,...,o) ^ nnd only if Ptv 
has a trivial kernel. 


11 




2. The time evolution is trivial on the 6-operators, i.e. 




( 2 . 22 ) 


Given this, by setting A{t) := Tt{A), understood component-wise as above, one finds 


A{t) = 


(2.23) 


3. The quadratic nature of S'at implies that any induced thermal state, i.e. p = j tr 
with /3 G (0, oo), or any eigenstate of Sn is quasi-free (cf. ||3|). The same applies to the 
functionals which result from these thr'ough decorations by the fermion parity operator, in 
particular, 

((•)) := tr {{■)Pnp) / tr {Pnp) (2.24) 

assuming tr (Pnp) A 0- 

The last observation will be essential in calculating the correlation functions (12.61) . 


2.2 Correlation functions 

Using the Jordan-Wigner transformation, the spin correlations functions (12.61) can be explicitly 
expressed in terms of correlation functions involving the a-operators. 


2.2.1 Reduction to Majorana correlations 

As a wai‘m-up, let us first consider all single spin correlations and all those involving Recall 
that we have set .A(f) = Tt{A) =: {af (f), o||"(f), • • • , o^(f), a)(f(f))^ the latter a notation we 
will use below. 

Lemma 2.1. Let p be a quasi-free state and assume [p, S'at] =0. One has that 

{a\) = {af) = 0, (a|) = i(a+ap , (2.25) 

and 


= (A(fr|)a2) = 0, 

- {a+{t)af) ■ (a^(f)a+), (2.26) 

for any 1 < < N and any f G M. 

Proof Inserting the Jordan-Wigner relation (12.71) we obtain for any 
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the final equality follows as these states are quasi-free. In fact, the expectation of the product is 
then a pfaffian, and there are an odd number of a-operators. The result for cr| is similar. The 
third identity in (12.251) immediately follows from = ia'^aj. 

For a derivation of (12.261) we proceed similarly using (12.71) and the fact that the state p is 
quasi-free: 






= i^ 



= 0 . 


(2.28) 


since the number of a-operators is odd. The result for cr^ is again argued similarly. 
To evaluate the remaining correlation, observe that 


('rt(o-l)o-i)) - (f)o+a^) -f (o+a^ ) • (a+a^) 

= (aj - {a^it)a-) ■ (o^(f)a+) (2.29) 


where the last equality again follows from the fact that p is quasi-free. In this case, the four-point 
function (o^ (f)a^(f)a+a“) reduces to a simple pfaffian which can be evaluated e.g. according 
to Wick’s rule (11.3b . Moreover, by time invariance we have (a^ (t)a^(f)) = (a^ aj). □ 

The cj^-correlation (12.26b is readily seen to decay in the distance |^ — ??| whenever the two- 
point functions involving the a-operators are known to do so. To establish a similar result for 
the correlations in the 12-plane, we again start from the Jordan-Wigner transformation (12.7b and 
write for w,w ^ { 1 , 2 }: 




'r}-l 


<') = n *4 apit)at-' J] ) ’ (^. 30 ) 


£=1 


where we introduced the abbreviation: 


\m=l 


#W := 


-|- w = \ 

- w = 2 


(2.31) 


The above average does not quite fit our needs if one aims to apply Theorem 11.31 or 11.41 We 
therefore rewrite the product 


= (det O) • ia^a^ ... Pn , 

using the identity (12.18b for the fermion parity. This brings the twisted average (12.24b into the 
equation: 

=(- 1 )-+-' i^+^-\detO) (trP^p) (-l)-'-i (2.32) 

X ((a^ {t)a^ (t) • • • a+_^(f)a^_^(f) a|^™(f)a^™'a++ia”+i ... a+a)^)) , 

for all 1 < ^ < 7/ < N, any f S M, and w,w' £ {1,2}. Here we have set b^/ = —#w'- In 
case (12.24b defines a quasi-free functional, the last expression (((...))) is the pfaffian 
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where x = (1,1, 2, 2,... , ^ — 1, ^ — 1, r/, r/ + 1, 7 / + 1,..., A^, N), the corresponding vector 

of signs is (+, bu,/, , +, —), and n = N + ^ — r]. This is precisely the 

setting of Theorem 1 1.3 1 or [T4l with distance given by 

r(x) = max \x 2 j -1 - X 2 j \ = - v\ ■ (2-34) 

j 

In order to apply these theorems, it remains to determine two-point functions associated with the 
a-operators. 

2.2.2 Calculating Majorana correlations 

All the relevant information concerning the spin correlations of interest is encoded in the fol¬ 
lowing 2N X 2N matrix: 

r^(t, s) := {A{t)A{sf) , t,seR. (2.35) 

If p commutes with Sn, then this correlation matrix only depends on the time difference, 
s) = — s, 0). Using (12.231) . it is clear that 

0 ) = A{0)A{0f) = 0)0^ (2.36) 

where we have similarly set r®(f, s) := We need only determine the static b- 

correlations and this is the content of the following 

Lemma 2.2. Assume either 

Case 1: p = j with /3 > 0, or 

Case 2: H]\[ has simple spectrum and p = |'I'q)('I'o| with a G {0,1}^, or 
Case 3: has a trivial kernel and p = tr with /3 > 0. 

Then for any t,s G M.' 

r^{t, s) = fp{HN), (2.37) 

where /p : M —)• M A the function given by: 

Case 1: fp{\) = 2(1 -y 
Case 2: fp{X) = 2xa„(A) 

with xAq characteristic function onto the set = {Xfaj = 0} U {—Xj\aj = 1}. 

Case 3: fp{X) = 2(1 - 

Proof Given (12.361) and the fact that Kn = AO, the claim is equivalent to showing that 

{BB'^)p = fpiiA ). (2.38) 
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One easily checks that, in the cases considered above, the off-diagonal expectations are zero, 
i.e., {h^h\)p = 0 for j / k, and so 


= »!)). = ® 


i=i 


i=i 


l-i{bjb. ) 


0 i 
J 'f’K-i 0 


(2.39) 


Here the last line results from explicit matrix multiplication using {bjY = 1 and b- bj = 
—b^bj. It thus remains to calculate i{b^b~)p in the cases mentioned above. 

1. In case of a thermal state, p = j tr we have 


i tr b 


j^-0Sn\ _ 


E 


N 


„-/3A,(2n,-l) 


{2nj 1 ) 


,-/3Afc(2n;,-l) 


v^nje{0,l} 

= — tanh(/3Aj) tre“^'^^ . 
0 i 


\nfc6{0,l} 


(2.40) 


This implies —i{b'^bj )p oj ~ tanh(i/3Aj). 

2. In case of an eigenstate, p = we have 

Yb'^bj)p = {^a, - l)^a) = 2aj - 1 


(2.41) 


Hence i{b'\b- )p{^ ■ 1) = {‘^cnj — l)sgn(iAj). Here that we require Xj / 0 for all 

j G {!,...,A^}, which due to the symmetry of the spectrum of Tf at is implied by the 
simplicity of the eigenvalues. From this the claim follows by distinguishing the cases 
aj G {0,1}. Note that the fact that has simple spectrum implies that Xa{H]\f) is 
well-defined. 

3. In case p = / tr , the calculation proceeds similarly to the first case. 

Since tr PNe~^^^ = (—2)^ sinh(/3Aj), we again need the assumption that Hn 

has a non-trivial kernel. 


□ 

Before turning to our main result, let us conclude this section with some historical remarks. 
Lemma 12.11 and 12.21 together with (12.301) or (12.321) yields a general expression for the time- 
dependent spin correlation functions of the XY model in terms of pfaffians involving the single¬ 
particle Hamiltonian Hn entering 

{afit)a!’pis)) = (<5f , (2.42) 

where {5|^} denotes the canonical orthonormal basis in ^^({1,..., A^}; C^). The fact that spin 
correlations in the XY -chain are expressible in terms of such pfaffians (or determinants) is an 
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observation which dates back to the seminal paper ifT/Tl for the time-independent case. In the 
homogeneous case these explicit expressions are used to show that the ground-state 

correlations in the 12 -direction exhibit an algebraic fall-off - a fact which should be conhasted 
to the exponential decay [230l below in the presence of an additional random field {t'g}. 

Explicit expressions for the time-dependent correlation functions go back to ll22]| (for the 3- 
direction) and ifT^ |71 (for the 12-direction). They have been the starting point for numerous 
further studies (see, e.g. 123). 


2.3 Main result: dynamical localization 

Our main result in this section concerns the case in which the one-particle Hamiltonian Hn 
is random and can be proven to exhibit shong-dynamical localization (cf. (11.181) '). A standard 
example of a random version of Tf at is the case that the spin coupling parameters {/r^} and { 7 ^} 
are constant and the external magnetic field {v^} forms iid random variables. We will discuss 
fhe applicabilify of fhe following general fheorem in fhis case below. 

Theorem 2.3 (Sfrong dynamical localization in spin chain). Suppose that the single-particle 
Hamiltonian H at associated with the spin-chain Sn is a random operator on £^{{1,, N }; ) 

which for all £ N satisfies: 

1. Hn has almost-surely simple spectrum. 

2. the eigenfunction correlator of exhibits complete strong dynamical localization in the 
sense that for a// G {1,..., N}: 


sup E 
#,befit} 


sup 

/6L°°(R) 

L||/l|oo<l 


{5*,f{Hn)5'’^) 




(2.43) 


with some N independent constants C,p, ^ (0, 00 ). 

Then the time-dependent spin correlations associated to either of the states corresponding to 
(1231) also exhibit strong dynamical localization in the sense that there are some C', p! G ( 0 , 00 ) 
for which, given any N G N.' 


E 


sup 

.teK 






< C' 


(2.44) 


for all w, w' G {1, 2,3} and all r/ G {1,..., A^}. 

Proof. In case re = 3 or re' = 3, fhe claim immediafely follows from Lemma lTTl and (12.421) wifh 
fp G bounded by ||/p||oo = 2, cf. Lemma |T2] In fhis case, fhe only non-hivial correlation 
is 


E 

sup (Tt(o-|)a^) - 

< 2 sup E 

sup 



LteM J 

#,befit} 

.teK 

- 


< 4(7e"^l^"'^l . (2.45) 
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In case w, w' € {1, 2}, we first restrict the discussion to the case of eigenstates, p = |'I'q)('I'q,| 
and envoke the representation (I2.32I) . Given (12.211) . the prefactor in (12.321) is bounded by one. 
One thus has for w, w' G {1, 2}: 




)CF„ 


{a"[{t)a^ (t) • • ap {t)arr '• ..a%aj^) 


(2.46) 

Since eigenstates are quasi-free, the right-hand side is the pfaffian (12.331) (with ((•)) replaced by 
(•)q,). The claim thus follows from Theorem [T3] using 


E 


sup sup 
a ieR 






(2.47) 


for all ^,7], Note that (12.471) follows from Lemma IT2] and assumption (12.431) . Moreover, as 
indicated in (12.341) . the distance of the configuration of Majorana fermions entering the pfaf¬ 
fian (12.331) is 1^ — /i|. 

In case of thermal states, p = j tre“^'^^, the result in case w,w' G {1,2} follows 

from the above, since 




< sup 




(2.48) 


The claimed bound is hence a consequence of Theorem 11.41 with the help of (12.471) and taking 
the first remark below Theorem 1 1.4 1 into account. □ 


Several remarks apply: 

1. As was shown in lH Prop. A.l], Sn and hence H]\] has simple spectrum for Lebesgue- 
almost all {u^} G M^. Taking {u^} independently distributed random variables with a 
single-site distribution which is absolutely continuous hence implies that H has almost 
surely simple spectrum. (Since the latter is symmetric about the origin, this in particular 
implies that the kernel of 77^^ is trivial almost surely.) 

2. As was explained in Subsection 12.11 in the isotropic case (7^ = 0) and for homogeneous 
spin coupling = p), the Hamiltonian 77jv reduces to (two copies of) the Anderson 
model with random potential In this case, strong dynamical localization in the sense 
of (12.431) is known to occur for iid random variables under fairly general conditions on the 
single-site distribution (cf. Q and references therein). 

In the non-isotropic, but homogeneous case (7^ = 7 and p^ = p) a. less complete picture 
is generally available. A result in ifT^ covers the regime of large disorder in case {u^} 
are iid with absolutely continuous distribution with a compact support. In ifTOl strong 
dynamical localization (12.431) is established for strong enough spin coupling \p\. 

Theorem 12.31 applies to all eigenstates and the thermal states. If one just aims at a localiza¬ 
tion statement concerning thermal states or ground-state, without any dynamics, less has to be 
assumed. 
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Theorem 2.4 (Localization for thermal states or the ground-state). Suppose that the single¬ 
particle Hamiltonian H]\f associated with the spin-chain Sn is a random operator on £‘^{{1,, Ai}; C^) 
which for all G N satisfies: 

1. Hn has almost-surely a trivial kernel. 


2. the Green function ofH]\i at zero exhibits fractional moment localization in the sense that 
for some s G (0,1) and all ^,r] € {1,..., A^}.' 


max sup E 
#,be{±} 


{SfdHN 





(2.49) 


with some N independent constants C,p, ^ (0, oo). 

Then the thermal spin correlations exhibit localization in the sense that there is C', p,' G (0, oo) 
for which, given any N G N.' 


E 






(2.50) 


for all w, w' G {1, 2,3}, all ry G {1,..., A^}, and all fi G (0, oo] . 


Proof. We will only give a proof in case G (0, oo) since the ground-state case (3 = oo follows 
by a limiting argument. 

In case m = 3 or u>' = 3, we proceed as in the proof of Theorem 12.31 In particular, in 
the only non-trivial case w = w' = 3, we use the first estimate in (12.451) in which /p(A) = 
2(1 -|- Using an argument from 121, we may write for Cjs = ^ G [1,1 -|- n/(3) and 

all A G C\{^ I n G Zodd} with | Im A| < e^: 


fp{X)=2Qi{X) + 2Q2{X) with 


Qi(A) := 


1 

¥ 


E 


nezodd 

r I n 

W 


/ zvrn 


-1 


TTlnl . 
20 


Q 2 {X) := / fp{u) 


u — iep — X u -\- icjs — X 


du . 


(2.51) 


(2.52) 


The contribution {5^^,Q2{Hn)5^^) is estimated with the help of a Combes-Thomas bound, 
cf. (11.161) and 13. The remaining contribution is estimated using (12.491 ). More explicitly: 


2E 


{6f,Qi{HN)d\ 


< 




E 

nezodd 

2/3 ^^P 


m 


l — S 




E 






n2-sf-i 1 




n=l 


(2.53) 
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where the last sum is bounded by a constant times (3^ for all /3 > vr. 
In the other cases w,w' G {1,2} we rewrite using (12.321) : 


=|trPAre ^^^1 


(2.54) 


X 


((a| 






an a, 


' 




• • • «Ar«Ar 


))/^ 


The arising pfaffian (12.331) satisfies the requirements of Therem ll.4l with Mq ^ = | tr Pjqe \ ^ 
0. To verify the other assumption (11.301) in this theorem, we note that by Lemma |2^ 

{{a*a\))p = (5f , 9p{\) = 2(1 - • (2.55) 

Since gp{X) = fp{—^ — A) we may use (12.511) together with the fact that 7rj{2f3) < ep and that 
the kernel of Tf at is trivial to rewrite 


qYHn) = 2Qi + 2Q2 - ^iv) . (2.56) 

The contribution of the second terms is again bounded using the Combes-Thomas estimate 
from 13. For its application note that ep — 7rj{2/3) > 1/2. The first term is estimated simi¬ 
larly as above. □ 


Some remarks: 


1. Quite generally, it is known that thermal states associated to one-dimensional, many-body 
quantum lattice systems satisfy exponential decay of correlations, or exponential clus¬ 
tering. In fact, Araki showed HI that analyticity arguments allow one to use Ruelle’s 
classical transfer matrix methods, see e.g. m, to prove that the Gibbs state of (e.g. finite 
range) one-dimensional systems satisfy exponential clustering at any positive temperature. 
Consequently, Araki’s result yields, deterministically, exponentially decaying bounds on 
thermal states of the XY-model. By contrast, our averaged bounds, in this random setting, 
are not only more explicit, they are also uniform, in the sense that they survive the /3 —)• oo 
limit. 

2. Some previous results concerning decay of correlations in random XY-models exist. In 
HU, some bounds on static correlations of certain observables in the ground state of the 
isotropic XY-model are considered. In lfT4ll . a bound on averaged, static ground state 
correlations is proven, again for the random isotropic XY-model. More precisely, for a 
chain of length > 1, a bound of the form 

]E(|(A5) - {A){B)\) < CN\\A\\\\B\\e-f^\^-^\ for all A £ B € (2.57) 

is obtained by combining a zero-velocity Lieb-Robinson bound, a Lifshitz tails estimate 
for the Anderson model, and well-known methods, see e.g. |[20ll . for deriving correlation 
decay in the ground state of gapped many-body systems. (Suggestions for improvement 
of this method can be found in |[T3l .l 
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3 Upper Bounds on Certain Bordered Determinants and 
Pfaffians 

The technical core of this paper, on which the proofs of our main results rest, are two estimates 
on certain bordered determinants and pfaffians. 


3.1 Determinants 


The following is the main new technical result for determinants. 

Theorem 3.1. Consider a complex matrix M £ ([;;;(«+i)x(”^+i) the following block structure 


j 

f “ 

T 

vi 

V2\ 


M = 

Wl 

A 

B 

(3.1) 

1 

[W2 

C 

D 



with a £ C, column vectors vi,wi £ and V 2 ,W 2 £ C'^, and blocks A £ B £ 

C £ and D £ with p + q = n. If \\M\\ < 1, then we have that 

|detM| < |a| + ||r; 2 || + \\wi\\ + ||i?|| + 2y^||ni||(||r(;i|| + ||i?||) (3.2) 


While this bound is not sharp, it also does not result from a straightforward application of 
Hadamard’s inequality |[T9l which asserts that 

n+1 

|detM| < TT \\rj{M)\\ < HMf min ||rj(M)|| . (3.3) 

j=t 

Here rj (M) denote the row vectors of the matrix and the last inequality follows since the Eu¬ 
clidean norm of any row is bounded by the matrix norm. Since we deal with matrices satisfying 
||M|| < 1, the determinant is then bounded by miuj ||rj(M)||. None of the row vectors, in 
general, have a norm comparable with the right side of (13.21) . In addition to Hadamard’s in¬ 
equality (13.31) . the proof of Theorem l3.1l is based on a change of basis and the invariance of the 
determinant under row operations. 


Proof of Theorem [01 Let (7 £ be a unitary transformation for which 


(C/^;i)^ = (0,...,0,||r;i||) (3.4) 

and take V £ to be a unitary which transforms AU'^ into an upper triangular matrix: 

/ai ... 


VAU^ = 


with some Uj £ C. Lifting these matrices 


(3.5) 


V 


OLr) 


U : = 


U 0 
0 


£(^nxn^ V: = 


V 0 
0 


£ C" 


(3.6) 
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we have 


M 


1 0 
0 V 


M 


1 0 
0 


a 

{Uvif 

\ 


Vwi 

VAU^ 

VB . 

(3.7) 

W2 

CU^ 

D 



Since the unitary transformations leave the norm as well as the modulus of the determinant 
invariant, an application of Hadamard’s inequality (13.31) yields 

I det M| = I det M| < min ||rj(M)|| . (3.8) 

3 

We now distinguish between two cases with a variational parameter e > 0. 

Case 1: Suppose that la^l < e. The norm of the (1 + p)th row of M can then be estimated as 

||rp+i(M)||= ^\{Vwi)p\^ + |ap|2 + \\rp{VBW < Ihill + e + ll^ll (3-9) 


where we used the fact that ||rj(yi?)|| < < ||f?||- Thus we have |detM| < 

||rui|| + e + ||i?||. 

Case 2: Suppose that Icepl > e. We then use row operations (which leave the determinant of M 
invariant) to eliminate the non-zero entry in Uvi from the first row. The other components 
in the first row are then modified as follows: 


a' ■.= a — s {Vwi)^ 

V 2 1-^ v '2 := V 2 — srp{VB ), wifh s := 


Pi 


(Xr\ 


(3.10) 


As a consequence, the norm of this modified first row can be estimated by 

+ lln^ll^ < la'I + Hn^H < |a| + |s| ||i^ri’i|| + 1111211 + l^l ||i’p(^-B)|| 


1 IPi II 

< |a| + - ||iii||||mi|| + ||i;2|| + (3.11) 


Therefore 


det M 


nil 


< min ||rj(M)|| < |a| + ||n2|| + ^ {\\wi\\ + \\B 


(3.12) 


Summarizing, the sum of the left sides of (13.91) and (13.121) constitute an upper bound on 
I detM|. Optimizing over e > 0, i.e., taking e = Y^||ni||(||t(;i|| + ||i?||), we then arrive at 
the bound claimed in (13.21) . □ 


21 












3.2 Pfaffians 


Since pfaffians belong to a less popular branch of linear algebra, let us start this subsection by 
reviewing some basic facts which will be of relevance. (For proofs and much more, see ifT^ .i 
Let M e ^ skew-Symmetric matrix, i.e. = —M. In the even case, i.e. m = 2n 

for some n G N, the pfaffian is defined by 

1 

^ ^ ^ sgn(7r) n a^( 2 j_i),^( 2 i) (3.13) 

Tr£S2n j = l 

where S2n is the symmetric group of permutations and sgn(7r) is the sign of the permutation 
71" £ <S' 2 n- (Taking the skew-symmetry into account this definition is seen to coincide with (ll.3l) .i 
The pfaffian of any skew-symmetric matrix with m odd is defined to be 0. It is also convention 
to define the pfaffian of a 0 x 0 matrix to be 1. 

Pfaffians share many similarities with determinants. First, they are invariant under certain 
elementary row operations which must be partnered with corresponding column operations to 
preserve skew-symmetry: 

1. Let M be the matrix obtained from M by multiplying a given row and the corresponding 
column of M by a constant A. Then pf[M] = Apf [M]. 

2. Let M be the matrix obtained from M by simultaneously interchanging two distinct rows 
and the corresponding columns. Then pf[M] = —pf[M]. 

3. Let M be the matrix obtained from M by taking a multiple of a given row and the cor¬ 
responding column and adding it to another row and the corresponding column. Then 
pf[M] = pf[M]. 

Next, pfaffians satisfy a Laplace expansion. The simplest case is an expansion along the first 
row/column, 

2m 

pf[M] = , (3.14) 

1=2 

where is the sub-matrix obtained from M by simultaneously removing two rows and two 
columns; namely those corresponding to 1 and i. 


Our new estimate concerns pfaffians of skew-symmetric matrices M G with 

the following block structure: 


M = 


/o 


V 


a V 


0 


1 
T 
wi 

A 


u 


TX 


T 

wi 

B 

c) 


(3.15) 


where a G C, the columns vi,wi G while V 2 ,W 2 G with p + q = n, and the blocks 
A G C2px2p, B G €2^x29, and C G C 29 x 2 g with both A and C also skew-symmetric. As the 
remainder of the matrix is determined through skew-symmetry, we leave it blank. 
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We will assume that M models a correlation matrix which, in particular, entails that the mod¬ 
ulus of its pfaffian is bounded. More generally, the following notion is tailored for our purposes. 


Definition. A skew-symmetric matrix M G ^‘inx 2 n ^ correlation structure 

of depth k G {0,1,... ,n} with constant Mq G (0, oo) if the pfaffian of all sub-matrices 
Afji which result from simultaneously eliminating the rows and columns labeled ji, j 2 ,... 

satisfy 


pfM: 


J1 J2i 


< Mo 


(3.16) 


for all disjoint integers ji,j 2 , ■ ■ ■ ,j 2 i G {1,..., 2n} and alH G {0,... , k}. The case ( = 0 by 
definition corresponds to no eliminations, i.e. the bound | pf M\ < Mq. 


We then have the following result. 

Theorem 3.2. Let M G be a skew-symmetric matrix which has the block- 

structure (13.151) and a correlation structure of depth 2 with constant Mq. Then 


2 p 


IpfMl < Mo j |a| -f ||i;2||i -f ||fi||i||i + ||fi||i||t(^2||i ^ \\rj{B)\\i 

i=i 


(3.17) 


where rj{B) G are the row vectors of B G and || • ||i denotes the 1-norm. 

The proof is based on two lemmas. The first is a straightforwai'd implication of the multi¬ 
linearity of pfaffians as expressed in (13.141) . 

Lemma 3.3. In the situation of Theorem \3.2\ 


|pf M| < (|a|-I-||u2||i-I-||ui||i||r(;i||i)Mo-I-||ui||i||ui2||i sup pf 

i6{l,...,2p} 
A;e{2p-|-l,...,2(p-|-i})} 


A B 
C 


j,k 

(3.18) 


Proof An application of the Laplace expansion (13.141) yields 

2p 2q 

pfM = api + ■ (3-19) 


f=l 


f=l 


Using the assumed correlation structure, the first and third term above are bounded by |a|Mo 
and ||u 2 II iMq, respectively. For the remaining sum, we again Laplace expand along the first row 


of"!,( 2 ^ 7 )^ 



/o . 


T—1 

II 


A B 

c) 


(3.20) 


(i+i) 


2q 


, 2 ,( 2 -ri),( 2 -rfe) 




k=l 


A B 




J21 
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with suitable exponents ej{k) e {0,1}. The first sum is bounded by ||rui||iMo. This gives rise 
to the third term in (I3.18I) . The second sum in is bounded by ||ru 2 ||i times the supremum in the 
right side of (13.181) . □ 

It remains to estimate the pfaffians which appear on the right side of (13.181) . 

Lemma 3.4. Consider a skew-symmetric matrix D G (hgfofjfi 


D = 


A B 
C 


(3.21) 


with blocks A G C(2p-i)x(2p-i), B G c(2p-i)x(2-?-i)^ and C G c(2'?-i)x(2<?-i). (Here all 
integers n,p,q > 1 and n = p-\- q — 1.) If in addition, D has a correlation structure of depth 1 
with constant K, then 

2p-l 

\pfD\<KY,\\rjm\i (3-22) 

j=t 

where rj{B) denotes the j-th row of B. 

Proof Since A G is skew-symmetric, its kernel is non-trivial. By the rank- 

nullity theorem, the range of A has dimension smaller or equal to 2p — 2, and consequently the 
columns of A and the rows of A are linearly dependent. Thus, there are //i,..., p 2 p-i S C, not 
all zero, for which 

2p-l 

^= 0 (3.23) 

where we have denoted by rj{A) the j-th row of A. As the pj do not all vanish, choose jo G 
{1,..., 2p — 1} satisfying \pj\ < \pjQ\ / 0 for all j. Without loss of generality, we will assume 
that jo = 1. Since the pfaffian is invariant under joint row/column operations, we may use (13.231) 
to eliminate the first row/column of D. In fact, 

/O 0^ 6f\ 

pfL> = pf( A Bj (3.24) 


where A = G 2 )x( 2 p 2 ) is the sub-matrix of A obtained by deleting the first row and 
column and B G (2ij 1) jg matrix with row vectors r 2 {B), ... , r 2 p-i{B) G 

Moreover, by inspection, it is clear that the first row is given by 

bl = ri{B) + ^^r,iB). (3.25) 

j>2 

We now Laplace expand along the first row on the right side of (13.241) and obtain 

2g-l 

0 . 26 ) 

k=l 
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Using the correlation structure for D and the fact that \fij\ < |/ii| we arrive at 

2p-l 

\piD\<K ll^illi <kY, \\rjm\i , (3-27) 

i=i 

which concludes the proof. □ 

Combining these two lemmas, there is a short proof of Theorem l3.2l 

Proof of Theorem \T2[ Using Lemma it is clear that we need only estimate the supremum 
on the right side of (13.181) . Each of the corresponding pfaffians is of a skew-symmetric matrix 
satisfying the assumptions of Lemma [3^ with K = Mq. □ 


4 Proof of decay of multipoint correlation functionals 

4.1 Proof of Theorems [HI and II .2! 

We organize the proofs of Theorems 11.11 and 1 1 .21 similarly. In each, for any n > 1, we consider 
a pair of configurations x = {xi,X 2 , • • • , Xn) and y = (yi, 2 / 2 , •'' > Un) both in ZP which we 
assume to be fermionic and ordered, cf. (11.121) . The configuation’s distance is attained at an 
optimizing pair jo G n}, i.e. 

D{x,y) = max \xj - yj] = \xj^ - . (4.1) 

l<j<n 

Without loss of generality, we will assume that 


Xjo < Vjo ( 4 - 2 ) 

since the roles of x and y may be interchanged in the case that Xj^ > yj^ . For convenience of 
notation we will relabel the particles 

jo !-)► 1, (1,..., jo - 1) 1-^ (2,..., jo), and (jo + 1,... , n) (jo + 1,..., n). (4.3) 


After this relabeling the n x n correlation matrix which we are interested in is given by M = 
with io{j, k) := (4^., p{sj,tk)5yj. It has the following structure: 


where we have set 


a := a;(l, 1) 


and 


A B 
C D 


fa Vi 

V2\ 

= 1 tui A 

B 

\w2 C 

d) 

/a;(l,2)\ 

/w(2,l) 

v={ : 

, m = : 


\w(n, 1) 

(w(j,A:))2<j.^ . 


(4.4) 


(4.5) 
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The sub-decomposition of the vectors v = (^ 1 ,^ 2 )^ and w = {wi,W 2 Y' and thereby also the 
matrix M into the blocks A, S, C, D is done according to the following rule in which we use 
the above relabeling 

a;(l, /c) is a component of vi if and only if ^ d{xi,yi ), 

oj{k, 1) is a component of wi if and only if ^ d{xi,yi). (4.6) 

This renders A and D square matrices. Moreover, we have the following estimates: 

Lemma 4.1. Let x and y be fermonic, ordered configurations in lA, i.e. such that holds, 
and set 5 := D(x, y). With respect to the relabeling introduced in d4.il) - d4.6l) above, one has 


V 2 \\ < p{xi,xi+tj, llwill < ^/5(yi - i,yi) (4.7) 

e>5/2 £>s 

- X] ^ y(xi - xi +/), (4-8) 

£>1 £'>5/2 


where p{xj,yk) := sup^ ^gj \{6xj, p{s,t)6yf^)\ and£,£' € Z. 

Proof. Recall that wlog (5 = yi — xi, see (14.21) . By the labeling rule (14.61) . we have that: 


1. The components of V 2 correspond to y^ > xi + 6/2. Thus, 


V2 


l(4i,y(si,4)5yfe)P < 

k: 

yk>xi+S /2 


p{xi,Vk)- 

k: 

yk>x\+&l2 


(4.9) 


Since the configuration y is assumed to be fermionic, i.e., yj 7 ^ y^ for all y / k, the right 
side above is trivially estimated by the first term in (14.71) . 


2. The components of wi correspond to yy < xi + 6/2 < yi. In this case, it must be that 
xy < xi. As a result, 

^ p{xpyi) (4.10) 

3- 

Xj<X\ 

similar to before. This is clearly bounded by the right side of (14.71) . 


3. The components of B are of the form uj{j, k) corresponding to y^ > xi + (5/2 and Xj < 
xi. Since the operator norm of B is trivially bounded by the Frobenius norm ||i ?||2 and 
moreover. 


Bh 


Y \{dxj,p{Sptk)6yf}\'^ < 

\ j: k: 

\j yk>xi+Sl2 


Y Y 

yk>xi+Sl2 


(4.11) 


(14.81) readily follows. 


□ 
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Using the above lemma and our general estimate for determinants (Theorem l3.ll) . we can now 
easily prove Theorems 11.11 and 1 1 .21 

Proof of Theorem \Tl\ Fix a pair of fermonic, ordered configurations x and y in IP. With re¬ 
spect to the relabeling described in (14.11) - (14.61) above, denote by M the corresponding correlation 
matrix which has the structure of (14.41) . Since ||/?(s,f)|| < 1, it is clear that ||M|| < 1. Using 
Theorem 13. 11 we have the estimate 

|detM|< \a\ + ||u 2 || + ||w^i|| + ||.B|| + 2y/||r’i||(||rui|| -f \\B\\) . (4.12) 

With Lemma |4~T] and the assumed decay, i.e. (11.141) . we have with 6 := D{x, y), 


a + ^2 + \\wi\ 


+ pll < +67^1 


-^.LK{l+5/2) 


e=o 


-f (7 ^ -h C ^ 


,-ij.K{e+e'+5/2) 


£=0 


£,£'=0 

oo 


oo \ 

< ^ ^ X] -f ^(1 -h I 

V (=0 1=0 ) 

< 4(7/(/io) . (4.13) 

Similarly, since ||ui|| < ||ri(M)|| < ||M|| < 1, the bound 


2Vlbi||(||w^i|| + II-BII) < 2a/||u;i|| -f ||B|| 

< 4V(7/(/io) 

follows. This completes the proof. □ 

Proceeding similarly, one has the following. 

Proof of Theorem 1772] We again denote by M = Mg^t the correlation matrix (making the time- 
dependence explicit). Arguing as above, the bound 


E sup IdetM^^tl <E sup -f E sup ||(us,t) 2 || +E sup ||(ms,t)i| 

\s,te/" / \s,tei I \sei,tei^ I \tei,sGi” 


-f E sup ||.Bs, 4 || -f 2, 
\s,tei” J \ 


E sup 


('l«s,t)l||+ sup ||S^,f| 
s,tei^ 


(4.15) 


readily follows from an application of Jensen’s inequality. 


With Lemma 14.11 and the a priori 
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estimate (11.191 ). we find that 


( \ i 6 ^ "2 S. 

E sup |detM,,t| < Ce-^^^ + C- - - + C- 




This completes the proof. 


1 — et* 1 — 






+ 2\ C 


o — ^lS 


+ C- 




(1 - e V l-e {l-e 


< 


SmaxIC, vCj 

- 1 —!- Lp A* 4 

(1 — e~>^Y 


(4.16) 

□ 


4.2 Proof of Theorem [lT3] and [TT4l 

Without loss of generality we will assume that the configuration x G is ordered according 
to (11.261) . There is some jo G {1,... ,n} such that r{x) = \x 2 j„ — X 2 j^-i\- We may now 
relabel jo ^ 1, (1 ,..., jo - 1) ( 2 , • • • , jo), and (jo + 1,. .., n) (jo + 1,..., re), such 

that the skew-symmetric 2re x 2re matrix featuring in Theorem 1 1.31 which has entries uj{j, k) := 

, can be assumed to have the following block structure 


M = 


^0 a uf V 2 \ 
0 W2 

A B 

\ c/ 


(4.17) 


with a := oj{l, 2 ) and 



( ‘^(1,3) \ 


^cu(l,2jo + l)\ 

/ cu(2,3) \ 

/a;(2,2jo + l)\ 

Vl = 


, V2 = 

: w'l = 


, u;2 = : 


yw(l,2jo)y 


y a;(l,2re) / 

ycj(2, 2jo)J 

\ a;(2,2re) y 


and 


^ ~ (^( 7 ) ^))3<i<fc<2n ■ 


(4.18) 


This decomposition has the property that the following norms are small (e t^^(''(*))) under the 
assumed decay of the two-point function: 


2n 2jo 

k=2jo+l k=3 

2jo 2n 

|77||cxd,oo 'y ^ y ^ ■ 

fc=3 l=2jo+l 


More precisely, we have 


(4.19) 

(4.20) 
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Lemma 4.2. For any Majorana configuration labelled such that (11.261) holds and r(x) = \x 2 j„ — 
X 2 j^-i\, we have 


OO 

e=o 

OO 

I'Wllll < 2j2Kx2ja,X2jo-l -i), 
e=o 


OO 

\V2\\1 <2j2pix2jo-l,X2j,+i) , 
i=Q 

OO 

|t(;2||i < 2 E Kx2jo,X2jo+^) 

e=o 


|^||oo,oo < 4 ^ Pix2jo-l -^,X2jo +^0 

e,e=o 


(4.21) 


where p{xj,yk) := max#^^g{±} | {5 %, p{s, t)5l^)\. 


Proof. Since there are no more than two Majorana Fermions (in fact of opposite flavor on 
each lattice site, the claim immediately follows from the definition of the vectors. □ 


Using the estimate (11.141) on p we may hence conclude (in a similar fashion as in (14.131) ): 


max{||ni||i, ||r(; 2 ||i} <2C'^e < 2CI{po ), 

£=o 

CX) 

max{||r; 2 ||i,||n;i||i} < < 2(U/(/ro) , 

OD 

Halloo,OO <^C ^-uK{r{x)+l+l') < /^CI{pq) e-^t^-t^o)K{r{x)) _ 

£,£'=0 

These estimates allows to apply Theorem l3.2l to the pfaffian of (14.171) and hence give a proof of 
Theorem 1 1.3 1 

Proof of Theorem [Q] As discussed above the skew-symmetric 2n x 2n matrix featuring in 
Theorem 1 1.3 1 has the block structure (14.171) . 


pf 


l<j<k<2n 


|pfM| . 


(4.23) 


Since (cc, r) is assumed to be a quasi-free pair and tc is bounded by Mq as a functional, the 
matrix M has a correlation structure of arbitrary depth with constant Mq, i.e. 


pfM: 


Jl J2v J2i 



# 


it. 


Jl/ 


.. .a. 


^21 


itj2l) 


^2 



< Mo . (4.24) 


Here we used the fact that eliminating rows and colums in the Pfaffian is equivalent to inserting 
the corresponding Majorana operators since {at{f))‘^ = 1. Moreover, the Majorana operators 
are bounded by one. We may hence apply Theorem 13.21 
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The norms of the bordering vectors have been estimated in (I4.22I) . Since the sums of the norms 
of the row-vectors of B are bounded according to ||rj(S)||i < ||-B||oo,oo> we conclude 


|pf M| < Mo + 2CI{fio) (1 + 2CI{fio)) 

. (4.25) 

This yields the claim. □ 


Proceeding similarly, one has the following. 

Proof of Theorem [7~?1 We again denote by M = the correlation matrix (making the time- 
dependence explicit). Arguing as above, the bound 


E sup — , ^ 
\ s,te/" -^0 


< E sup Mg ® I pf Mg^tl 3 


< E sup \as,t\ +E sup ||(us,t) 2 ||i 




1 
3 


-fE sup ||(us,t)i||i E sup ||(u;s,f)i||i 


. se/.ie/" 




-fE sup ||(us,t)i||i E sup ||('u;s,t) 2 ||i E sup \\Bs,t\ 


(4.26) 


. se/,4e/" 


.sei,tei^ 


, s,tei^ 


follows from the bound | pf Mg^tl < -^o and Holder’s and/or Jensen’s inequality. The expecta¬ 
tion values in the right side may now be estimated using Lemma 1431 and the assumption (11.281) . 
Note that since the maximum in (11.281) is over a finite set, we may assume that this bound also 
holds for the maximum inside the expectation value. □ 
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